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Testing the stability of fundamental constants 
with the 199 Hg + single-ion optical clock 

S. Bizefj S.A. Diddams, U. Tanaka, C.E. Tanner,Q W.H. Oskay, R.E. Drullinger, T.E. Parker, 
T.P. Heavner, S.R. Jefferts, L. Hollberg, W.M. Itano, D.J. Wineland, and J.C. Bergquistf] 

Time and Frequency Division, National Institute of Standards and Technology, 325 Broadway, Boulder, CO 80305, USA 

(Dated: February 2, 2008) 

Over a two-year duration, we have compared the frequency of the 199 Hg + 5d 10 6s 2 Si/2(F = 

0) < > 5d 9 6s 2 2 D 5 / 2 (F = 2) electric-quadrupole transition at 282 nm with the frequency of the 

ground-state hyperfine splitting in neutral 133 Cs. These measurements show that any fractional 
time variation of the ratio vcs/vn s between the two frequencies is smaller than ±7 x 10 _15 yr~ 1 
(lcr uncertainty). According to recent atomic structure calculations, this sets an upper limit to a 
possible fractional time variation of gcs(rn e / 'm p )a 6 ' at the same level. 

PACS numbers: 06.30.Ft, 32.30.Jc, 32.80.Pj 
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The development of string theory models aiming at 
a unified description of gravity and quantum mechanics 
has renewed interest for improved experimental tests of 
Einstein's Equivalence Principle (EEP). Indeed, a com- 
mon feature of these models is that they allow for, or 
even predict, violations of EEP Q. These include vi- 
olation of the universality of free-fall as well as varia- 
tion of fundamental constants with time and space. In- 
terestingly, a recent analysis of the spectrum of quasars 
H suggests that the fine-structure constant a may have 
changed over the cosmological timescale (10 10 yr). The 
Oklo reactor analysis j| || on the other hand puts a 
stringent limit to possible variation of a on the geolog- 
ical timescale (10 9 yr). Owing to their high accuracy, 
comparisons between atomic frequency standards based 
on different atomic species and/or types of transitions 
provide one of the best ways to perform laboratory tests 
of the stability of fundamental constants. Present and 
future efforts to improve atomic frequency standards, in 
both the optical and the microwave domains, will im- 
prove these tests, leading to significant constraints on 
theoretical work aimed at a unified theory. 

In this letter we describe frequency comparisons con- 
ducted over a two-year period between a 199 Hg + single- 
ion optical clock and a 133 Cs fountain atomic clock that 
set a new stringent limit to a possible variation of funda- 
mental constants. The theoretical background for such a 
test is given first. We then describe the experiment and 
conclude with the results of the test. 



the frequency i^Hg 

„2 2] 



of the 



1!)!) 



Hg 4 



In the experiment, 
5d 10 6s 2 S 1/2 (F = 0) < — > 5d 9 6s 2 2 D 5/2 (F = 2,m F = 0) 
electric-quadrupole transition at A = 282 nm is compared 
to the frequency vcs of the ground-state hyperfine tran- 
sition 6S 1/2 (F = 3, m F = 0) < — ► 6S 1/2 (F = 4, m F = 0) 
in neutral 133 Cs. Including relativistic and many-body 
effects, j/jjg can be expressed as i^Hg — Ry Fn g (a), where 
Ry = RoaC is the Rydberg constant expressed as a fre- 
quency, and Fn g (a) is a dimensionless function of the 



fine-structure constant a = e 2 /Aire o he. Similarly, the 
hyperfine frequency of cesium can be approximated by 
"Cs — gcs(rn e /m p )a 2 R y F Cs (a), where g Cs is the gyro- 
magnetic ratio of the 133 Cs nucleus and m e /m p the 
electron-to-proton mass ratio. Fn s (a) and Fc s (a) are 
calculated in @, 0|. With numerical values from Jq], we 
find 



a ^-lnF Hg (a) ~ -3.2 
oa 

a ^-lnF Cs (a) -+0.8. 
oa 



(1) 
(2) 



Therefore, sequential measurements of the ratio vcs/^Rg 
actually test the stability of the product of fundamental 
constants 5cs(w e /m p )a 6 . Note that the sensitivity to 
variation of the individual constants needs to be known 
to only 1-10 % in order to adequately describe how this 
comparison constrains possible variations of these funda- 
mental constants. It is therefore justified to omit higher- 
order terms in the expression of the hyperfine splitting 
of 133 Cs, such as finite nuclear size. The high sensitiv- 
ity of the Hg + vs Cs frequency comparison to a change 
of a arises from the large relativistic effects encountered 
in heavy atoms combined with the negative sign of the 
relativistic effects in 199 Hg + . These factors make the 
i99jjg+ 2 S\/2 < — ► 2 £>5/2 optical transition one of the 
best choices for a search of a variation of a [0] . 

The 199 Hg + single-ion optical frequency standard has 
been described previously [M || and only the main fea- 
tures are outlined here. A single 199 Hg + ion is stored 
in a radio-frequency (rf) spherical Paul trap held at a 
cryogenic temperature (4.2 K). It is laser-cooled to near 
the 1.7 mK Doppler limit using the strongly allowed 
2 S 1/2 (F = 1) < — ► 2 P 1/2 {F = 0) transition at 194 nm, 
as shown in Fig. |l|. After a cooling and state prepara- 
tion phase (~ 16 ms), the ion is left in the 2 Si/ 2 (F = 0) 
lower state of the clock transition. The cooling laser is 
then switched off and the probe laser light at 282 nm is 
directed onto the ion for a typical duration of T = 50 ms. 
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FIG. 1: Partial level scheme of 199 Hg + . The 194 nm 

2 Si/ 2 (F = 1) < ► 2 Pi/2{F = 0) transition is used 

for Doppler cooling, state preparation and detection. The 
282 nm electric-quadrupole transition from the ground state 
2 Si/ 2 (F = 0) to the metastable 2 D 5/2 (F = 2,m F = 0) state 
provides the reference for the optical clock frequency. 



The cooling laser is turned on again to determine the 
ion's internal state using the technique of electron shelv- 
ing |], The 194 nm fluorescence photons emitted by 
the ion are counted for 15 ms. The absence of scattered 
photons indicates that the ion has been excited to the 
2 D5/ 2 (F = 2) state by the 282 nm probe laser. Simi- 
larly, a typical count rate of ~ 6000 s^ 1 from the ion 
indicates that it has remained in the 2 Si/ 2 (F = 0) state 
after the probe period. When the ion is detected in the 
2 D$/2 excited state, the 194 nm is left on until scattered 
photons are detected again, indicating that the ion has 
spontaneously decayed to the 2 S\/ 2 state. A new inter- 
rogation cycle is then started. 

The 282 nm radiation used to probe the clock transi- 
tion is obtained by frequency-doubling 563 nm light from 
a dye laser in a deuterated ammonium dihydrogen phos- 
phate (AD*P) crystal, as shown in Fig. |^. The light 
from the dye laser is stabilized to a resonance of a sta- 
ble high-finesse (T ~ 200 000) Fabry-Perot cavity. Two 
acousto-optic modulators (AOMl,AOM2) driven by pre- 
cision rf sources shift and fine-tune the frequency of the 
563 nm radiation so as to match the frequency of the 
282 nm radiation to the S — D resonance frequency. By 
compensating the predictable linear drift of the stable 
reference cavity (~ 500 mHz s _1 ) using AOM1, we real- 
ize an interrogation oscillator with a fractional frequency 
instability of 3 x 10~ 16 between 1 s and 10 s, corre- 
sponding to a linewidth of 640 mHz at 282 nm Q. For 
a typical probe time T = 50 ms, which is shorter than 
the 90 ms spontaneous decay time of the 2 D 5 / 2 state, 
the observed linewidth is Fourier-transform limited with 
a full width at half maximum of ~ 16 Hz, much larger 
than the frequency fluctuations of the probe laser during 
the probe pulse time. Under these conditions, the fluctu- 
ations of the measured transition probability, and hence 
the frequency stability of the Hg + frequency standard, 
are limited mainly by atomic quantum projection noise 



In our trap, the radial and axial secular frequencies are 
1.1 MHz and 2.2 MHz respectively. Both frequencies are 
much higher than the recoil frequency v = h/2m\ 2 ~ 
12.6 kHz, where m is the ion mass and h is the Planck 
constant. The observed spectrum therefore consists of 
a central feature at the 2 S 1 / 2 (F = 0) < — > 2 D 5 / 2 (F — 
2,m,F = 0) transition resonance frequency, which is free 
of the recoil frequency shift, together with vibrational 
sidebands at integer multiples of the trap secular fre- 
quencies ]l0|] . The spectrum also exhibits features corre- 
sponding to the carrier and vibrational sidebands of the 
2 S 1/2 {F = 0) <— ► 2 D 5/2 (F = 2,m F = ±1,±2) "Zee- 
man components." Typically, the frequency difference 
between the clock transition and the first Zeeman com- 
ponents is ~ 6 MHz, which corresponds to a quantization 
magnetic field B w 0.3 mT. 

In order to lock the frequency of the probe laser to 
the atomic resonance, a square-wave frequency modula- 
tion is applied to the 563 nm beam using AOM2 so as 
to probe the carrier of the clock transition alternately 
on each side of the resonance. Typically, 24 measure- 
ments are averaged on each side of the resonance to find 
the transition probability. The difference between these 
two transition probabilities reveals the detuning between 
the center frequency of the probe laser and the center 
of the atomic resonance. This error signal is used in a 
digital servo loop that steers the average frequency of 
the probe laser to the center of the atomic resonance, 
applying frequency corrections to the synthesizer driv- 
ing AOM1. The time constant of this servo loop is on 
the order of Ti 00 p <~ 15 s. The analysis of the frequency 
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FIG. 2: Experimental setup for the absolute measurement 
of the frequency of the 199 Hg + optical clock in terms of the 
SI second defined by the ground-state hyperfine splitting of 
133 Cs. 
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corrections indicates that for r > ri oop , the fractional 
frequency instability of the probe light stabilized to the 
atomic resonance is 5 — 7 x 10 -15 t^ 1 / 2 , in good agree- 
ment with the quantum-limited instability expected for 
our experimental conditions (trap secular frequencies, ion 
temperature and measurement cycle time). 

In order to compare the Hg + optical frequency stan- 
dard to other frequency standards, some fraction of the 
563 nm light at the output of AOM1 is launched into 
a 180 m long optical fiber and delivered into a separate 
room where it is frequency shifted by AOM3. AOM3 is 
used to actively cancel the optical path length fluctua- 
tions of the fiber in order to preserve the high degree of 
coherence of the light [|l2[ |l3) . 

The frequency of the 563 nm light is measured with 
an optical frequency comb generated by a mode-locked 
TkSapphire laser whose femtosecond pulses are spec- 
trally broadened in a microstructure fiber 0, |l5[ [l7| . 
The optical spectrum at the output of the fiber consists 
of equally spaced, phase-coherent modes with frequen- 
cies /„ = nf r + fo, where f r is the repetition rate of 
the mode-locked laser, n is an integer and fo is a fre- 
quency offset caused by the difference between the phase 
and group velocities in the laser cavity. In our setup, 
the repetition rate is typically f r ~ 1 GHz and the spec- 
trum spans more than one octave from ~ 520 nm to 
1170 nm. The repetition rate is detected directly with a 
fast photodiode. The offset frequency fo is detected by 
the self-referencing method, where the frequency-doubled 
red part of the comb 2/„ = 2nf r + 2f is heterodyned 
with the blue part of the comb / 2n = 2nf r + f . We 
also detect the beatnote ft between the 563 nm light 
from the Hg + optical frequency standard and the closest 
mode of the optical frequency comb f m = mf r + fo (m 
is known from previous coarse measurements of the Hg + 
S — D transition frequency). In practice, fo and fb are 
phase-locked to precision rf sources by acting on the cav- 
ity length and the pump power of the femtosecond laser, 
respectively. Finally, f r is measured by counting the low- 
frequency beatnote between f r and a third rf synthesizer. 

In order to perform absolute optical frequency mea- 
surements in terms of the SI second, all synthesizers and 
frequency counters are referenced to the 5 MHz output 
of a hydrogen maser with a typical frequency instability 
~2x 10 -13 at 1 s. The maser itself is part of an ensem- 
ble of 5 masers and 3 commercial cesium clocks used to 
realize the local NIST atomic time scale, which is in turn 
periodically calibrated using the NIST-F1 cesium foun- 
tain primary standard [ Q , as well as international ce- 
sium standards. The frequency of the reference hydrogen 
maser is known within 1.8 parts in 10 15 with respect to 
the ground-state hyperfinc splitting of 133 Cs. As shown 
in recent investigations jl9|, , the additional noise and 
inaccuracy from the optical frequency comb itself is neg- 
ligible at this level. Using this setup, it is thus possible 
to perform absolute measurements of the frequency of 



the Hg + optical standard (together with the frequency 
of each component of the optical comb) with a fractional 
frequency uncertainty of 1.8 x 10~ 15 . Typically, measure- 
ments are performed for 2 hours, leading to a statistical 
(type A) fractional frequency uncertainty of 2.4 parts in 
10 15 , which corresponds to a 2.5 Hz uncertainty on the 
frequency of the 282 nm stabilized light. 

The full evaluation of all systematic effects of the Hg + 
optical standard is still under way. At the present time, 
10 Hz is a conservative upper bound for the total sys- 
tematic (type B) uncertainty, as shown by the following 
preliminary analysis of systematic shifts. The second- 
order Zeeman frequency shift of the clock transition is 
given by 5v z = K^B 2 , Kf = -18.925(28) kHz mT" 2 
|pi[ . The typical field B w 0.3 mT corresponds to a 
second-order Zeeman frequency shift 8fz ~ —1.7 kHz. 
Therefore, the 1.5 x 10 -3 fractional uncertainty on co- 

(2) 

efficient K z leads to a 2.6 Hz uncertainty on the clock 
frequency. For the same value of the bias magnetic field, 
the sensitivity of the clock frequency to field fluctuations 
is 11 Hz fiT . In our unshielded environment, i?-field 
fluctuations up to ±0.2 fiT have been observed, leading 
to a 2.2 Hz uncertainty. The most troublesome frequency 
shift is the electric-quadrupole shift 8vq due to the cou- 
pling between the atomic electric-quadrupole moment in 
the upper 2 D 5 / 2 (F — 2, nip — 0) clock state with electric 
field gradients due, for example, to stray charges on the 
trap electrodes. From our trap geometry and the bias 
voltages applied to compensate for stray electric fields, it 
is possible to calculate an upper bound for an electric field 
gradient at the location of the ion and to estimate that 
8vq < 1 Hz [^l|. In addition, two successive versions 
of the trap electrodes (gold and molybdenum surfaces) 
have been used, leading to no detectable change of the 
clock frequency. Similarly, the background pressure of 
helium (the only species remaining with significant pres- 
sure at 4.2 K) has been changed by more than one or- 
der of magnitude (estimated by measuring the ion heat- 
ing rate due to background collisions) without producing 
any detectable frequency shift. Empirically, this implies 
that the He background pressure shift is smaller than 
1 Hz. At a temperature of 300 K, the black-body radia- 
tion shift is —0.08 Hz, and it is considerably lower in the 
4.2 K cryogenic environment [|TJ. At the Doppler cool- 
ing limit of 1.7 mK, the second-order Doppler shift due to 
thermal motion is —0.003 Hz. Finally, by suitably com- 
pensating the stray static electric field, systematic shifts 
related to the trapping oscillating electric field (second- 
order Doppler shift due to micromotion, AC Stark shift) 
are made smaller than 0.1 Hz p2f . 

Over a period of two years, we have performed 20 
measurements of the 199 Hg+ S — D transition frequency 
with respect to the 133 Cs ground state hyperfine split- 
ting. Figure [| shows these measurements, corrected for 
the second-order Zeeman frequency shift, with their sta- 
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FIG. 3: Absolute frequency measurements of the 1 Hg 4 



S U2 {F = 0) 



Db/2(F = 2) transition with respect 



to the 133 Cs ground state hyperfine splitting defining the SI 
second. The plot shows the deviation of each measurement 
from the weighted average value with its statistical ±1<t er- 
ror bar. The horizontal axis is the Modified Julian Date 
(MJD 52582 =Nov. 4, 2002). 



tistical (type A) error bars. The absolute frequency of 
the Hg + optical standard is given by the weighted av- 
erage of these data: m g = 1 064 721 609 899 143.4 Hz. 
The total statistical uncertainty is only 1.0 Hz. Measure- 
ments in Fig. |^ clearly show a reproducibility better than 
10 Hz at 1.06 x 10 15 Hz, the most stringent comparison of 
optical and microwave frequencies to date. The total un- 
certainty of the measurement is dominated by the 10 Hz 
systematic uncertainty of the Hg + clock. With a 10~ 14 
fractional uncertainty and a measurement period of two 
years, our measurement constrains a possible fractional 
variation of ^Cs/z-'Hg at the level of ±7 x 10 -15 yr _1 (la 
uncertainty) . 

This result can be interpreted as constraining a pos- 
sible fractional variation of gcs(w e /mp)a 6 ' at the same 
level: ±7 x 10~ 15 yr . Assuming that any change in 
this quantity is due to the a 6 factor, we derive an upper 
bound for a possible linear variation of the fine-structure 
constant: \a/a\ < 1.2 x 10 -15 yr -1 , a factor of 30 im- 
provement over j^] . However, there may be a signifi- 
cant change of gc s or m e /m p due to possible variation of 
the strength of the strong and electroweak interactions. 
In fact, recent theoretical work argues that within the 
framework of a grand unified theory a fractional varia- 
tion of a is necessarily accompanied by a fractional vari- 
ation of m e /m p that is ~ 38 times larger p3. Fortu- 
nately, comparisons between two optical clocks will test 
the stability of a alone M. One interesting possibility is 
to compare the 199 Hg + S — D transition to the 657 nm 
1 So{m = 0) < — > 3 Pi (to — 0) transition in neutral 40 Ca 
fli] , . With an independent constraint to the stability 
of a, comparisons involving hyperfine transitions, such 



as the present fn g (optical) vs ^cs comparison or the 
z/Rb (microwave) vs ^c s comparison |25|, ^6| , will test the 
stability of the strong and electroweak interactions. Ref- 
erence [ p7| investigates in more detail possible laboratory 
tests of the stability of fundamental constants. 
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